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ABSTRACT 

Context. The Taurus-Auriga association is perhaps the most famous prototype of a low-mass star forming region, surveyed at almost 
all wavelengths. Unfortunately, like several other young clusters/associations, this T association lacks an extensive abundance analysis 
determination. 

Aims. We present a high-resolution spectroscopic study of seven low-mass members of Taurus- Auriga, including both weak-lined and 
classical T Tauri stars designed to help robustly determine their metallicity. 

Methods. After correcting for spectral veiling, we performed equivalent width and spectral synthesis analyses using the GAIA set of 
model atmospheres and the 2002 version of the code MOOG. 

Results. We find a solar metallicity, obtaining a mean value of [Fe/H]=-0.01±0.05. The Q--element Si and the Fe-peak one Ni 
confirm a solar composition. Our work shows that the dispersion among members is well within the observational errors at variance 
with previous claims. As in other star forming regions, no metal-rich members are found, reinforcing the idea that old planet-host 
stars form in the inner part of the Galactic disc and subsequently migrate. 

Key words, stars: abundances - Galaxy: open clusters and associations: individual: Taurus-Auriga 



1. Introduction 

The Taurus-Auriga (Tau-Aur) association is the nearest large 
star-forming region (SFR), with a close distance of d=14Q pc 
(e.g., Kenyon et al. 1994). Given its lack of massive O/B stars, 
this T association in the past 50 years has become a standard 
region to study the low-mass star formation processes. Tau-Aur 
is widely spread across the Northern sky (~100 square degrees, 
Scelsi et al. 2007a) and contains hundreds of low-mass mem- 
bers; Rebull et al. (2010) found 148 new candidate members, of 
which 34 were confirmed by spectroscopic follow-up. The bulk 
of stars presents an average age of ~1 Myr (e.g., Briceno et al. 
1999). 

The whole region has been observed at almost all wave- 
lengths, from infrared to X-rays (e.g., Itoh et al. 1996; Briceno 
et al. 1999; Davis et al. 2008; Luhman et al. 2010). Kenyon & 
Hartmann (1995, hereafter KH95) provided a thorough investi- 
gation of the stellar population, from class (proto-stars) to class 
II/III (i.e. classical and weak-lined T Tauri stars, respectively - 
hereafter CTTs and WTTs). They analysed optical and infrared 
photometric observations and presented colour-magnitude and 
colour-colour diagrams, luminosity and mass functions, along 
with information on the near-infrared excess and accretion prop- 
erties. Lithium abundance studies of this association have been 
also performed over the past decade, such as, e.g., Basri et al. 
(1991), Magazzu et al. (1992), Martin et al. (1994), and Sestito 
et al. (2008), who presented a (re )asse ssment of Li i abundances 
in WTTs and CTTs. Giidel et al. ('2007) carried out an extensive 
X-ray survey (XMM-A^ewfo« Survey of the Taurus Molecular 
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Cloud -XEST) covering an area of ~ 5 square degrees and con- 
centrating mainly on the higher stellar density regions. 

Whereas all these studies have helped to characterize the 
stellar population, Tau-Aur shares with other nearby SFRs the 
lack of an accurate abundance analysis. 

The determination of chemical composition of SFRs is in- 
stead critically important to a variety of astrophysical issues, in 
both planetary and stellar contexts, as we previously discussed 
in our pilot project focusing on the Orion association (D'Orazi 
et al. 2009, hereafter D09). At variance with OB associations, 
whose abundances can provide an independent test of the so- 
called triggered star formation scenario (e.g., Blaauw 1991) and 
indicate whether there is local enrichment, we would expect to 
measure a very homogeneous composition for T associations. 

As is well known, giant planets are preferentially found 
around metal-rich main sequence stars (Santos et al. 2004, 
and reference therein). The natural question hence arises as to 
whether in the first epochs of planetary formation (discs of T 
Tauri stars are commonly accepted as planet birthplaces) the 
probability of hosting a giant planet depends on the star's metal- 
licity. 

On the other hand, and very interestingly, several studies 
(e.g., Luhman 2004) have found significant differences between 
the initial mass function (IMF) derived for the Tau-Aur and those 
derived for denser systems containing massive members, e.g., 
the Trapezium, with the former containing a too small number 
of brown dwarfs according to the standard IMF (Briceno et al. 
2002) and a surplus of late-K and early-M stars (Luhman et al. 
2009). These differences can be attributed to different conditions 
of the environment where the stars were formed. Some studies 
have claimed that the IMF slope may depend on metallicity, in 
the sense that metal-rich environments tend to produce more 
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low-mass stars (M < O.IMq) than metal-poor systems (Larson 
2005; da Rio et al. 2009 and references therein). The obvious 
question is: could metallicity play a role in the lower fraction of 
very low-mass stars and brown dwarfs detected in Taurus with 
respect to Trapezium (Briceno et al. 2002)? More generally, 
we may ask whether there is a difference between chemical 
compositions of low-mass and high-mass SFRs or, whether, 
SFRs represent a chemically homogeneous class of objects? 

Available information on Tau-Aur abundances comes from 
the previous works by Padgett (1996) and Santos et al. (2008). 
First, Padgett (1996) analysed eight WTT members and derived 
a solar composition, i.e. an average value of [Fe/H]=0.01+0.13 
(rms), with a significant scatter among members ranging from 
-0.11 to -1-0.22. Santos et al. (2008) presented the iron abun- 
dances of three Taurus members: in this case, the spread in 
metallicity is also large, varying from -0.18 to +0.05 with a 
mean value of [Fe/H]=-0.07+0.12 (rms). Given the quite large 
star-to-star variation obtained by both studies, a new accurate 
abundance analysis for this association is warranted, to mini- 
mize the impact of the observational uncertainties and ascertain 
whether the metallicity scatter among members is indeed real or 
not. 

In a different framework, Scelsi et al. (2007b) obtained coro- 
nal abundances for a sample of X-ray bright members of Taurus, 
finding that the iron abundance in the corona is significantly (by 
a factor of ~ 5) lower than the solar photospheric value. More 
generally, their coronal abundances revealed a pattern in agree- 
ment with those obtained in previous studies. Specifically, X- 
ray observations of young and/or active stars have unveiled a 
so-called FIP (first ionization potential)-dependent abundance 
trend: low FIP elements (FIP<10 eV, like iron, silicon, nickel) 
are under-abundant in the corona with respect to the photo- 
spheric values (the so-called inverse-FIP effect, for further de- 
tails see Brinckman et al. 2001; Audard et al. 2003; Scelsi et al. 
2007b). As stressed by Scelsi et al., however, this FlP-related 
pattern for Tau-Aur, as well as for other young stars, was ob- 
tained by comparing the coronal abundances to the solar photo- 
spheric values, rather than the stellar photospheric abundances. 
Different or no FIP effects were indeed found when considering 
the correct stellar photospheric abundances (e.g., Sanz-Forcada 
et al. 2004; Maggio et al. 2007). Determining the photospheric 
abundances of Tau-Aur members with measured coronal abun- 
dances thus appears critical to confirm whether a FlP-related ef- 
fect is present in these stars or not. 

In this work, we present a high-resolution spectroscopic 
study based on seven low-mass members of the Tau-Aur asso- 
ciation. The sample was selected specifically to determine ac- 
curate abundances, by discarding, e.g., fast rotators, strong ac- 
cretors, and/or binary stars. Most importantly, this is one of the 
first sample that also includes T Tauri stars that have not yet dis- 
sipated their circumstellar discs and reveal moderate accretion 
from the surrounding material. Details of the observations and 
data reduction are given in Section |2l Section [3] describes the 
estimate of and correction for spectral veiling, along with the 
abundance analysis procedure. The results and the scientific im- 
plications are illustrated in Sections |4] and |5j while in Section|6] 
a conclusive summary is given. 



2. Observations and data reduction 

Our sample consists of seven low-mass members of the Tau-Aur 
association, of which two are CTTs and five WTT stars, respec- 



tively. The sample was selected from the large database provided 
in XEST by Gudel et al. (2007). 

Single K-type members were chosen with projected rota- 
tional velocities (v5/«i) not larger than ~30 km s"' and a mass 
accretion rate below log M=-7.5 Mq yr"'. For stars rotating 
faster than 30 km s"' the abundance analysis becomes infeasi- 
ble, because of rotational blending. The upper limit to the accre- 
tion rate allowed by including in the sample also CTTs and, at 
the same time, discarding extremely strong accretors, for which 
abundance analysis is no longer possible. Basic properties of our 
sample are given in Table [1] where we list the star name, the 
observation date, and the exposure times in Columns 1, 2, 3, re- 
spectively; signal-to-noise ratio (S/N) per pixel measured around 
the Li I line at 6708 A in the co-added spectra are reported in 
Column 4. The V magnitudes and spectral-types are listed in 
Columns 5 and 6, respectively, while the initial effective tem- 
perature (Teff) from Giidel et al. (2007) is given in Column 7. 
Finally, in Column 8 we show the classification as CTT (Class 
II) or WTT (Class III) stars, while projected rotational velocities 
from XEST catalogue are in Column 9. 

The observations were carried out in service mode at the 
Telescopio Nazionale Galileo (Roque de los Muchachos, Canary 
Islands) with the high-resolution spectrograph SARG, during 
seven nights of observations between August and December 
2007. We used the 61 A mm^' CD3 yellow grism (Aiaze=5890 
A), which resulted in a spectral coverage between 4620 A and 
7920 A , along with a 1 .6 slit, providing a resolution R ~29000. 
As we discuss below, we also acquired a solar spectrum with 
the same instrumental set-up, which permits us to carry out a 
differential abundance analysis (see Sect. O with respect to the 
Sun. Data reduction was performed within the echelle context in 
MIDAS following the standard procedures of bias subtraction, 
flat-field correction, order definition and extraction, sky subtrac- 
tion, and wavelength calibration. Spectra of the same star ac- 
quired in different exposures were then co-added, after checking 
that there were no variations in radial velocity. The spectra of 
seven of our sample stars are shown in Figure [1] in the spectral 
window between 6700 and 6760 A, where several Fe i features 
are present. 

3. Analysis 

3.1. Estimate of veiling and equivalent-widtli correction 

Spectral veiling, which often affects spectra of pre-main se- 
quence (PMS) stars, was estimated in the same fashion widely 
described in D09. Briefly, we compared the equivalent widths 
(EWs) of the sample stars with those of stars with similar pa- 
rameters belonging to the older clusters IC 2602 (~ 30 Myr) and 
IC 2391 (~ 50 Myr): the veiling parameter is expressed as r = 
(EWic,/EWtau)-1. 

Nine spectral lines, of different atomic species, were em- 
ployed in this computation, namely Ca i 5857.5 A, Ca i 6102.7 
A, Ca I 6122.2 A, Fe i 6546.3 A, V i 6624.8 A, Ni i 6643.6 A, 
Fe 1 6662.5 A, Ca 1 6717.7 A, and Ti 1 6743.2 A. In a subsequent 
step, we then computed for each star the mean r value coming 
from the different spectral lines; these average values, along with 
the standard deviation, are listed in Table[3] All the WTTs of our 
sample have negligible veiling, as expected; on the other hand, 
the two CTTs, GM Aur and BP Tau, exhibit quite large r val- 
ues. In particular, our values of 0.23(+0.09) and 0.60(±0.07) are 
in excellent agreement with the previous estimates by Basri et 
al. (1991), who found 0.25(±0.1) and 0.60(+0.1), respectively. 
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Fig. 1. Portion of the spectra of our seven sample stars between 6700 A and 6760 A; the dashed hne separates the WTTs from the 
two CTTs (GM Aur and BP Tau). 



Table 1. Log of observations and stellar properties. 



Star 


Obs.Date 


Exp.time 




V 


Sp.Type 


TefF 


Class 






(yy-mm-dd) 


(sec) 




(mag) 




(K) 




(kms-') 


GM Aur 


2007-09-08 


3600 


110 


12.03 


K7 


4060 


II 


12.4 


BP Tau 


2007-12-20 


3600 


110 


12.16 


K7 


4060 


II 


7.8 






1168 














V826 Tau A 


2007-09-19 


3600 


90 


12.11 


K7 


4060 


III 


4.2 


Hubble 4 


2007-09-20 


2x3000 


120 


12.67 


K7 


4060 


III 


12.8 






2150 














RXJ0438.2+2023 


2007-10-30 


3600 


110 


12.22 


K2 


4900 


III 


8.0 






2550 














RXJ0423.7+1537 


2007-09-19 


2200 


140 


11.29 


K2 


4900 


III 




LkCa 19 


2007-12-20 


1270 


100 


10.85 


KO 


5250 


III 


18.6 



for the two stars. Finally, the correction for spectral veiling was vate communication) Q. We performed both EW and spectral 
applied to the two "veiled" stars, by multiplying for (1 + r) the synthesis analyses, using the drivers abfind and synth, respec- 
observed EWs to obtain the correct values. lively. Radiative, Stark, and collisional broadenings were treated 

in the standard way, adopting for the last one the classical Unsold 

(1955) approximation. 

3.2. Abundance analysis 



LTE abundances were derived by applying the spectral analy- 
sis program MOOG of Sneden (1973, 2002 version) and using 
the GAIA set of model atmospheres (Brott & Hauschildt, pri- 



' 'http : //www .hs . uni- Hamburg . de/EN/For/ThA/p hoenix/index.html 

Biazzo et al. (2010) provide a wide discussion of why this set of model 
atmospheres were chosen for these young, cool stars. 
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Table 2. Solar abundances derived from our analysis for Fe, Si, 
and Ni. As a comparison, we also show in Column 3 the standard 
abundances by Anders & Grevesse (1989). 



Element 


This work 


AG(1989) 


log n(Si i) 


7.54±0.03 


7.55 


log n(Fe i) 


7.51 ±0.06 


7.52 


log n(Ni i) 


6.24+0.02 


6.25 



3.2.1 . Equivalent width analysis and stellar parameters 

For the five slowest rotators in our sample (vsin/<18 km s"'), 
we derived iron abundances by measuring the EWs. The line list 
was based on the one provided by Randich et al. (2006) to which 
we added four iron lines from the list by D09; we refer to those 
papers for further details on atomic parameters and, in general, 
for the whole procedure. 

We derived solar abundances by analysing a SARG spec- 
trum of the Sun taken with the same resolution and grating of 
our sample stars (see Table |2]for solar abundances of Fe, Si, and 
Ni). Our complete line list includes 48, 2, and 13 features of Fe i. 
Si I, and Ni i, respectively; however, depending on the star tem- 
perature and S/N, we were able to measure a different number of 
lines for the different stars. The final number of lines used, after 
the 2cr clipping, is given in Table [3] Assuming Tgff 0=5770 K, 
log gQ-4.44, and ^©=1.1 km s"' (see Randich et al. 2006), we 
derived log n(Fe 1)0=7.5 1 ±0.04. The EWs for all the stars were 
measured using the IRAF task SPLOT and performing both a 
Gaussian fit and a direct integration to the line profile. In the case 
of Hubble 4, the bluer specti-al region (4600 A -6200 A) could 
not be included, since this spectral window was missed because 
of technical problems during the observations; the analysis was 
thus performed using the upper (i.e., redder) orders including 
only 13 iron lines. Along with iron, we also derived the abun- 
dances of Si I and Ni i for the five slow rotators for which EW 
analysis could be carried out. 

Initial stellar parameters were retrieved from the XEST 
source catalogue. In particular, we adopted Teft values derived 
from spectral types and the calibration of KH95. After apply- 
ing a 2cr clipping to the initial line list, we derived the final 
Teff values by zeroing the slope between log n(Fe i) and exci- 
tation potential (x). To infer microturbulence (^), we first as- 
sumed 1.4 km s"' as an initial value, and then imposed the 
condition of there being no trends between log n(Fe) and EW 
strengths. Surface gravities were computed using the expression 
log ^=4.44 H-log (M/M0)-log (L/Lo)-H4xlog Teff-15.0447, tak- 
ing mass and luminosity values from summary tables by Giidel 
et al. (2007). We could not spectroscopically optimize these 
gravity values, because of the lack of suitable Fe 11 lines in the 
spectra of the cool stars. For the warmest star within the slow- 
rotating sample, i.e. RXJ0438. 2-1-2023, we were able to mea- 
sure only one Fe 11 line, lacking statistics. Therefore, photometric 
gravities are adopted as final values. 

In Table [3] we list the final stellar parameters: the initial ef- 
fective temperatures agree very well with the final ones, having 
a maximum difference of 70 K (for BP Tau). The final microtur- 
bulence values range from 1 .3 km s ' for the coolest star BP Tau 
to 1.9 km s"' for the warmest one, namely RXJ0438.2-H2023. 



3.2.2. Spectral synthesis 

For the two WTT stars with vi/«i larger than 18 km s"', EW 
measurements were infeasible because of strong line blending. 
We thus derived a metallicity estimate using spectral synthe- 
sis method in a wavelength window of ~ 15 A, from 6700 
A to 6715 A, using the line list already employed in D09. We 
first optimised that line list by changing the \oggf values, when 
necessary, to obtain the closest agreement between our SARG 
solar spectrum and the standard iron abundance by Anders & 
Grevesse (1989). To compare the different observed spectra for 
these rapidly rotating stars with the synthetic spectra, we con- 
volved the latter with both a Gaussian profile corresponding to 
our resolution of R=29000 (FWHM=0.23 A) and a rotational 
profile, taking into account the Umb-darkening coefficients. 

The following stellar parameters were adopted. Values of Tefr 
and log g were retrieved from XEST; we then assumed micro- 
turbulence of ^=1.6 km s ', since this is the mean value ob- 
tained from the EW analysis of other sample stars. Finally, spec- 
tral synthesis allowed us to infer, as a by-product, the projected 
rotational velocities vs/ni (Table|3]l. 

In Figure s|2] and [3] examples of spectral synthesis are shown 
for the stars RXJ0423.7+1537 and LkCal9, respectively. They 
are discussed in Section |4] 



3.3. Error budget 

Internal (random) errors, affecting the abundances derived by 
the EW analysis, are due to (i) uncertainties in the oscillator 
strengths and EW measurements, and (ii) uncertainties in the 
adopted set of stellar parameters. The first type of uncertainty, 
i.e., log gf, should be minimised, given the differential nature of 
our analysis; on the other hand, uncertainties in EWs are clearly 
represented by cr[ in Table[3] Errors due to uncertainties in Tgff, 
log g, and ^ were estimated as described in D09, to which we 
refer for a complete and detailed description of the approach 
followed. In brief, we assumed typical errors of 60 K in Tetf, 
0.2 dex in log g, and 0.2 km s"' in ^ since larger errors in Tgii 
and ^ would introduce trends in abundances with x and EW, re- 
spectively. For gravity, we instead adopted the same uncertainty 
as D09. We then derived the variation in the iron abundance and 
both [Ni/Fe] and [Si/Fe] by changing one stellar parameter at 
the time, while keeping the others unchanged, and finally we 
quadratic ally added the three contributions to obtain the total er- 
ror (0-2). Furthermore, for the two "veiled" stars, we computed 
the abundance error due to the veiling estimate (0-3) which im- 
plied that, by assuming a conservative error equal to the rms of 
the mean r value (see Column 2 in Table |3]l, the variation in the 
iron abundance is +0.09 and +0.06 for GM Aur and BP Tau, 
respectively. 

Uncertainties related to best-fit determination and stellar pa- 
rameters affect the abundances derived from spectral synthesis: 
errors in the best-fit are around 0. 1 dex (tri). Errors in stellar pa- 
rameters were computed as described by D09: a variation of +60 
K in TefF and +0.2 km s ' in ^, will result in abundance changes 
of +0.06 dex; the effect of log g is instead smaller, reaching a 
maximum of ~0.04 dex. 

Systematic (external) errors, caused for instance by the code 
and/or model atmosphere, should not influence largely our final 
abundances, as widely described in D09. 
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Fig. 2. Left panel: Best-fit synthetic spectrum ([Fe/H]=0.00) of the star RXJ0423. 7-1- 1537. Right panel: Zoom on the small spectral 
window containing two Fe i features. Synthetic spectra obtained with the best-fit [Fe/H] value and its + 1 cr margins of error 
superimposed on the observed spectrum for comparison (solid, dotted, and dashed lines, respectively). 




Fig. 3. Observed and synthetic spectra of the star LkCal9. The "formal" best-fit is shown in both panels as a thick solid line (see 
text -Section l4.2l for a discussion on this case). 



4. Results 

4.1. Mean abundances 

The final abundance values are reported in Table [3] for both 
the EW and spectral synthesis analyses. For each star (Column 
1), we list the veiling parameter r in Column 2, while the final 
(adopted) set of stellar parameters ( Teff, logg, ^) are presented 
in Columns 3, 4, 5. In Column 6, we show the number of lines 
used in EW analysis (upper panel) and the vsini (km s"') as a 
by-product of spectral synthesis (lower panel); finally, metallic- 
ity [Fe/H] is given in Column 7 along with its corresponding 
errors. 



Given the difference between the methodologies, we com- 
puted two different average values: the first technique (EW anal- 
ysis) inferred [Fe/H]=-0.02+0.06, while considering the faster 
rotator RXJ0423. 7-1-1537 (and not considering the less secure 
estimate for LkCa 19, see Section |42]l we obtain a mean metal- 
licity of [Fe/H]=-0.01+0.05. The small standard deviation from 
the mean, which is much smaller than expected from the formal 
errors, indicates a homogeneous solar iron abundance for Tau- 
Aur with no star-to-star scatter. 

The a-element Si i along with the Fe-peak one Ni i 
also confirm solar values (see Table |4]i: the average ratios 
[Si/Fe] =0.04+0.03 and [Ni/Fe] =-0.02+0.05 indicate that Tau- 
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Aur members are characterised by the same abundance pattern 
as thin disc stars of similar metalhcity. 

Our study provides the first abundance determination using 
the EWs of "deveiled" CTTs, following the work of Biazzo 
et al. (2010) for the Orion complex. A moderately accreting 
PMS star was also present in the Padgett's sample, namely CV 
Chamaeleonis; however, after the estimate of spectral veiling 
(r=0.2+0.3), the EWs correction was not performed and the au- 
thor concluded that the derived iron abundance should be treated 
as a lower limit. We find that CTTs and WTTs share the same 
chemical composition; while this result indicates that the pres- 
ence of a circumstellar disc does not affect the stellar photo- 
spheric abundances, it also strongly supports the validity of our 
method for estimating spectral veiling. 

4.2. Comparison with previous works 

Only one of our stars is in common with previous studies in the 
literature: LkCal9. For this star, we obtain a formal best-fit with 
[Fe/H]=-i-0.2±0.2 (see TableO, suggesting a possible iron over- 
abundance. However, the uncertainties in the stellar parameters 
of this star prevent us from definitively assessing its metal-rich 
nature. As Figure |3] shows, even the formal best-fit synthetic 
spectrum cannot simultaneously reproduce both the Fe i lines 
at 6703 and 6705 A, suggesting that the adopted Teff, as de- 
rived from spectral-type (Tgff =5250 K), is likely to be over- 
estimated. We indeed note that for LkCal9 Martin et al. (1994) 
derived a value of Tefr=4343 K, which is probably too low. An 
intermediate value might be the correct one and would yield a 
lower [Fe/H]. A lower effective temperature would also yield a 
lower gravity. Since the available spectrum, and in particular the 
presence of several blended features, does not allow us to inde- 
pendently derive stellar parameters, at the present time, no final 
statement can be drawn about the metallicity of this star. 

This star was also previously analysed by Padgett (1996), 
who obtained a Tetf=5214+146 K, an iron abundance 
[Fe/H]=-0.08+0.08, and a significantly higher microturbulence 
^=2.8+0.5 km s ' from an EW analysis. The Padgett's line list 
contains about one third of spectral features with EWs larger 
than 150 mA with a maximum value of 330 mA : these high 
values are probably due to blended features, given the high rota- 
tional velocity of ~20 km s In other words, an EW analysis is 
really uncertain in this case, when weak metallic lines blend to- 
gether and the resulting abundance probably has to be revised. To 
conclude, given the uncertainties in eflFective temperature, hence 
metallicity, the star LkCal9 was discarded in the computation of 
the average abundance for Tau-Aur, as reported in Section|4T| 

Focusing on the overall metallicity, in Figure |4] our [Fe/H] 
determination in Tau-Aur is compared with the previous esti- 
mates by Padgett (1996) and Santos et al. (2008); their averaged 
values of [Fe/H]=-0.01+0.13 and [Fe/H]=-0.07+0.12, respec- 
tively, are in good agreement with our new estimate. Moreover, 
the [Si/Fe] =0.04+0.04 and [Ni/Fe] =-0.07+0.02 ratios derived 
by Santos et al. (2008) are similar to the values inferred here. 

However, the figure clearly shows that the [Fe/H] distribu- 
tion from our measurements is much narrower than those ob- 
tained by both those studies. As mentioned in Sect. [1] in the 
work by Santos et al. (2008) the metallicity ranges from -0.18 
to +0.05, while in Padgett (1996) the internal dispersion is even 
larger, with [Fe/H] ranging from -0.16 to +0.22 dex. Padgett 
concluded that, beyond the major uncertainties affecting young 
stars with respect to the older ones, "the standard deviation of 
[Fe/H] within the T Tauri star association is larger than older 



clusters". However, our homogeneous abundance determination 
shows that, at least within our sample, the star-to-star scatter is 
smaller than the observational errors and that previously detected 
variations probably reflect these uncertainties. The same con- 
clusion has been drawn by D09 and Biazzo et al. (2010), with 
both studies focusing on the abundance determination across the 
Orion complex. On the other hand, we note that the number of 
SFRs chemically analysed by considering a large sample of stars 
is still very limited; further investigations are necessary to con- 
firm (or refuse) the presence of real internal scatter. 



5. Discussion 

5. 1 . Tlie metallicity of Taurus-Auriga 

The distribution of [Fe/H] for the few SFRs and young associa- 
tions with available abundance measurements was discussed by 
Santos et al. (2008) and Biazzo et al. (2010). None of them is 
more metal-rich than the Sun and, in general, they appear to be 
a part of a rather homogeneous sample, characterised by similar 
iron abundance patterns. Biazzo et al. obtained an average metal- 
hcity [Fe/H]= -0.06 ± 0.04, considering SFRs within 500 pc 
of the Sun and including the value for Tau-Aur derived in the 
present paper 

Nevertheless, we note that small (but larger than the errors) 
differences in the metallicity of the regions are present, with the 
Orion Nebula Cluster (ONC) and Chamaeleon being the most 
metal poor (with [Fe/H]=-0.11 + 0.08 and -0.11 ± 0.11, re- 
spectively) and Tau-Aur, with a solar metallicity, being the most 
metal rich. Although the current sample is still too small to per- 
form reliable statistics and draw definite conclusions, it seems 
that low-mass SFRs do not behave differently from high mass 
SFRs, that metallicity is not related to the characteristics and po- 
sition of the region, and that the shape of the IMF is not related 
to the metallicity. In particular, we notice that a difference in 
metallicity is most likely not the reason for the different shapes 
of the IMF in the Chamaeleon and Taurus; indeed, as mentioned 
above, Chamaeleon is more metal poor than Taurus and, accord- 
ingly, should produce a smaller number of low-mass stars. We 
observe instead a peak at ~ O.8M0 in the IMF of Taurus, with 
a deficit of lower mass stars, at variance with other regions, like 
Chamaeleon, showing an IMF peak at lower masses (around 0. 1 
Mq, see Luhman 2004 and references therein). 

5.2. Coronal abundances 

Scelsi et al. (2007b) determined the coronal abundances of about 
20 T Tauri stars in the Taurus Molecular Cloud. They found a 
clear inverse-FIP effect, with median Fe and Si abundances ~0.2 
times the solar photospheric abundances. However, as stressed 
by Scelsi et al. and in Sect. 1 of this paper, the determination of 
photospheric abundances of young (active) stars is critically im- 
portant to comparing in a homogeneous way these values with 
the coronal ones and drawing more definitive conclusions on 
possible FlP-related trends. 

Three of our stars (BP Tau, V826Tau, and Hubble4) are in 
common with the sample of Scelsi et al. For all three stars, we 
obtain solar Fe and Si abundances from our analysis. This con- 
firms a large depletion of iron and silicon in the hot coronae of 
these stars with respect to the photospheres, which may be at- 
tributed to condensation of iron in grains, as suggested by Scelsi 
et al. (2007b). 
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Teff 


logg 




Nr lines 


[Fe/H]±cri±o-2±cr3 








(K) 




(km s"') 




(EWs) 


GM Aur 


0.23H 


t0.09(C) 


4100 


3.9 


1.7 


17 


0.04±0.08±0.08±0.09 


BPTau 


0.60H 


t0.07(C) 


3900 


3.7 


1.3 


18 


-0.07+0.05+0.09+0.06 


VS9fi Tan A 






4.000 




1.4 


78 

Z.O 


oo+n 00+0 08 


Hubble 4 






4000 


3.3 


1.5 


13^ 


-0.08±0.05±0.08 


RXJ0438.2+2023 






4850 


4.5 


1.9 


27 


0.03±0.09±0.05 














\ sinv 


[re/rlJ+(T 1 +cr 2 














(km s"') 


(Synthesis) 


RXJ0423.7+1537 






4900 


4.3 


1.6 


25±2 


0.00±0.20±0.06 


AVERAGE 














-0.01±0.05 


LkCa 19 






5250 


4.2 


1.6 


21+3 


0.20±0.20±0.06 



" Only the upper CCD was available 
Table 4. [Si/Fe] and [Ni/Fe] ratios for a sub-sample of our stars 



slar 


LSi/l-cJ±(r|±(r2±(r3 


Nr lines 


LNi/l-cJ±(r|±(r2±(r3 


Nr lines 


GM Aur 


0.04±0.09±0.07±0.03 


1 


-0.05±0.04±0.03±0.02 


4 


BPTau 


0.08±0.07±0.08±0.02 


1 


0.06±0.07±0.04±0.02 


5 


V826Tau A 


0.02±0.02±0.07 


1 


0.01+0.01+0.03 


10 


Hubble4 






-0.02+0.14+0.03 


2 


RXJ0438.2+2023 


0.03±0.07±0.06 


2 


-0.06+0.04+0.03 


4 


AVERAGE 


0.04+0.03 




-0.02+0.05 





5.3. Metallicity and planet formation 

Surveys of old solar-type stars have clearly shown that the prob- 
abiUty of hosting giant planets increases with the metalUcity of 
the parent star: the frequency of giant planets around stars of 
twice solar metallicity is about ~ 30%, to be compared with 
~3% for stars characterised by solar or sub-solar iron content 
(e.g., Santos et al. 2004; Fischer and Valenti 2005). The origin 
of this correlation between metallicity and planet incidence is 
still debated, with the coexistence of different hypotheses, from 
detectability (i.e., the higher the metallicity, the greater the like- 
lihood of planet migration and hence planet detection), to envi- 
ronmental effects (i.e., metals enhances the formation of dust, 
planetesimals, and finally cores). We mention in passing that the 
correlation between planet and metallicity seems to have two 
major caveats: (1) the different behaviours of dwarfs and giants 
(see Pasquini et al. 2007), giant stars hosting planets not appear- 
ing to be more metal-rich on average than star without planets; 
and (2) this correlation is no longer valid for [Fe/H] values rang- 
ing from -0.7 dex to -0.3 dex (see Haywood 2009). 

The high metal abundances of planet-bearing stars seem 
to be primordial (e.g., Pinsonneault et al. 2001; Ecuvillon et 
al. 2006; Gilli et al. 2006). A natural question is therefore 
whether young PMS stars with super-solar metalUcity exist, 
since their circumstellar discs are commonly assumed to be the 
planet birthplace. As mentioned in Sect. 5.1, all the young clus- 
ters/associations for which an abundance analysis has been per- 
formed to date do not contain any high-metallicity star. For six 



SFRs, Santos et al. (2008) derived a slightly sub-solar metal- 
licity, concluding that metal-rich stars in the solar proximity 
(and especially planet-host stars) might have formed in the in- 
ner Galaxy and then suffered radial migration across the Galactic 
disc. This scenario also seems to be confirmed by several studies 
based on both abundances and kinematics, such as Ecuvillon et 
al. (2007) and Haywood (2008). 

Given the quite limited sample, no final conclusion can 
be drawn about the possible presence of metal-rich SFRs, or 
at least, a population of metal-rich T Tauri stars inside them. 
However, we can speculate that at the present time in the solar 
neighbourhood, for its young clusters/associations, the probabil- 
ity of forming giant planets is not as high as in the inner disc. 
This does not mean that no giant planets can form (and possi- 
bly survive) at all in these young stellar aggregates, but rather 
that the frequency of giant planets could be somewhat lower 
than in the inner regions of our Galaxy, where the planet-host 
stars are suggested to originate. To explain this, two hypothe- 
ses have been proposed: (1) The higher metal content of the in- 
ner Galaxy and thus the higher dust-to-gas ratio can enhance the 
probability of forming gas giant planets (for details of the core 
accretion model for giant planet formation see, e.g.. Pollack et 
al. 1996; Ida & Lin 2004; Mordasini et al. 2009a,b). (2) As an 
alternative, metalhcity-independent approach, Haywood (2009) 
suggested that giant planet formation could be favoured in re- 
gions where the density of the molecular hydrogen (which is the 
primary constituent of these objects) is significantly higher. We 
note that the molecular ring of our Galaxy is located precisely at 
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Our work 




-0.2 0.2 



[Fe/H] 



Fig. 4. Comparison of our [Fe/H] (filled histogram) with previous estimates by Padgett (1996, empty histogram) and Santos et al. 
(2008, dashed histogram). 



the Galactocentric radii of 3-5 kpc where these metal-rich stars 
are observed (Haywood 2009). 



6. Summary and conclusions 

We have presented elemental abundances of seven stars, both 
CTTs and WTTs, belonging to the Tau-Aur association. We have 
found a very homogeneous, solar metallicity for this T associa- 
tion, deriving a mean value of [Fe/H]=-0.01+0.05. Both [Si/Fe] 
and [Ni/Fe] also exhibit solar ratios and agree very well with the 
observed abundance pattern of thin disc stars at the same metal- 
licity. In contrast to the previous determination of abundances in 
Tau-Aur, which spanned a wide range in [Fe/H], we conclude 
that the internal dispersion in metallicity for this association is 
very similar to values derived in older open clusters. Along with 
our previous project in Orion and the available estimates in the 
literature, all the SFRs surveyed to date seem to share a simi- 
lar chemical composition, suggesting a uniform ISM in the solar 
surroundings at the present time. In this context, we note that 
no metal-rich members have been detected in all the analysed 
young associations: this could validate the idea that metal-rich 
planet-host stars in the solar circle were formed in the inner disc 
and subsequently moved to their current location. Their not-so- 
high metallicity might reflect the rather low frequency of giant 
planets in the solar neighbourhood, at variance with that of the 
inner Galactic regions. However, given that only a small number 
of SFRs have been chemically characterised so far, and a small 



number of stars per region, further investigations are mandatory 
to definitely help us address this controversial, but intriguing is- 
sue. 
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